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ABSTRACT 

Bifurcation  to  high-dimensional  hyperchaos  is  observed  in  a 
driven  coupled  pendulum-flexible  rod  system.  When  the  rod  is 
in  resonance  with  the  pendulum,  the  system  changes  from  a  low 
dimensional  attractor  to  a  high  dimensional  attractor  abruptly.  It 
is  shown  that  high  dimensional  chaotic  dynamics  is  hysteretic, 
and  exhibits  extreme  sensitivity  with  respect  to  small  parameter 
changes.  Such  sensitivity  poses  a  problem  in  obstructing 
predictability  in  mechanics.  A  brief  discussion  of  sustaining 
chaos  to  prevent  resonance  is  included. 

INTRODUCTION 

One  of  the  simplest  problems  studied  in  nonlinear  physics  is  that 
of  the  forced,  damped  pendulum  (Baker,  1996),  (Miles,  1993). 
Studied  extensively  in  isolation,  the  pendulum  is  always 
attached  to  some  support,  such  as  a  stiff  rod,  and  is  used  in  many 
active  control  mechanisms,  such  as  vibration  absorbers 
(Starretl,1995),  (Cuvalci,  1996).  If  the  rod  is  sufficiently  stiff, 
one  expects  the  pendulum  dynamics  to  be  slightly  perturbed 
from  the  ideal  infinitely  stiff  case.  We  review  some  new 
dynamical  behavior  in  the  numerical  simulations  of  a  forced, 
damped  pendulum  coupled  to  a  linear  rod  which  is  flexible.  Wc 
examine  the  system  when  it  is  operating  in  a  resonant  mode, 
where  the  pendulum  frequency  is  half  that  of  the  fundamental 
frequency  of  the  rod,  as  well  as  a  non-resonant  case.  It  is  known 
that  when  the  rod  is  sufficiently  stiff,  (lie  dynamics  resides  on  a 
global  slow  invariant  manifold;  i.c.,  the  rod  is  slaved  to  the 
motion  of  (lie  pendulum  (Georgiou  and  Schwartz,  1999).  As  a 
result  of  this  slaving  motion,  the  dynamics  is  a  perturbation  of  a 
parametrically  driven  pendulum. 

However,  when  operating  at  resonance,  there  exists  some 


critical  amplitude  of  the  driving  force  that  causes  a 
discontinuous  change  from  periodic  behavior  to  hyper-chaotic 
behavior,  where  there  are  two  or  more  positive  Lyapunov 
exponents. 

In  contrast,  non-resonant  dynamics  exhibits  multiple  chaouc 
attractors,  one  which  is  statistically  constrained  to  a  low 
dimensional  surface,  and  one  which  is  noL  In  both  cases, 
multiple  attractors  as  well  as  long  lived  chaotic  transients  exist. 
Such  multiplicity  gives  rise  to  an  obstruction  in  predicting  the 
type  of  asymptotic  behavior  due  to  small  uncertainties  in 
parameters  and/or  data  (Schwartz,  ct  al,  1 999).  Either  chaotic  or 
periodic  behavior  may  result  for  given  initial  data.  To  eliminate 
the  possibility  of  periodic  resonant  behavior,  the  results  of  a 
new  parametric  control  procedure  to  sustain  chaos  arc  presented 
(Schwartz  and  Triandaf.  1 996),  (Triandaf  and  Schwartz,  2000). 

MODEL  EQUATIONS 

The  dynamical  system  we  consider  (see  Fig.l)  consists  of  a 
planar  pendulum  of  length  Lp ,  mass  M  ,  and  viscous 
dissipation  coefficient  Cp  attached  to  the^lower  end  B  of  a  linear 
viscoelastic  rod  of  length  Lr ,  cross-section  Ar ,  mass  density  pr , 
elasticity  modulus  Er ,  coefficient  of  internal  viscous  dissipation 
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Figure  1.  Rod-Pendulum  Geometry. 

The  upper  end  A  of  the  rod  is  subject  to  a  periodic  motion  x.(t). 
Let  u(x,  t)  and  6(r)  denote  respectively  the  displacement  field  of 
e  rod  with  respect  to  A  and  the  angular  displacement  of  the 
pendulum  mass. 

Jf5  =  f’T  =2V  derivative  with  respect  to  xis  denoted  by 
a  dot  ,  where  a>p  _  £  is  the  uncoupled  pendulum  natural 
equency  denote  normalized  spatial  and  temporal  variables. 
Letting  U  =  —  denote  the  dimensionless  displacement  field  for 
the  rod,  the  (dynamics  of  the  rod-pendulum  configuration  are 
described  by: 

6  =  -[l-VB’(x)  +  JP4(x)jsin0_25  6 

2  2 

^  4  V,TT^’  =  T)  +  2C/ltFT^C4.  X)  -  (1) 

2  2 


SPATIO-TEMPORAL  BEHAVIOR 
Resonance 


— ' - -  ^Iiuauun  wuere  the  rod  has  a 

fundamental  natural  frequency  of  2,  and  the  drive  frequency  i 
also  equal  to  2.  (Recall  the  natural  frequency  of  the  penduluL 
unity,  making  F  =  ± .)  Hie  other  fixed  parameters  of  the 
normalized  model  are  the  mass  ratio  (unity),  and  damping 
coefficients.  The  variable  control  parameter  is  forcing 
amplitude  a  We  plot  in  Fig.  2  a  bifurcation  diagram  of  the 

fofthe'I  ^  eXpami°n  t0  ^  (J-2)  is  observed 
for  the  angular  vdoeity  0f  the  pendulum  sampled  at  the  period 

ilt^ 


a 


Figure  2.  Bifurcation  diagram  of  pendulum  velocity 
sampled  at  the  forcing  frequency  as  a  function  of  a  .The 
coupling  parameter  is  p  =  ~  . 


V(0,x)  =  o 
2  2 

V^d.x)  =  -  7'(0,  0,  x)  ] cos 0 

.2 

7(6,  0,  t)  =  0  +  [1  -  A^(x)]cos0 


(2) 


The  details  of  the  model  and  parameter  scalings  are  given  in 
(Georgiou  and  Schwartz,  1999).  The  physical  parameters  that 
control  die  coupling  jjftween  the  pendulum  and  the  rod  are  the 
frequency  ratio  M  =  ^  and  the  mass  ratio  p  .  For  fixed  p  die 
lima  of  coupled  systeifi  Eqs.  (1-2)  as  g  ^  0  describes  the 
motions  of  the  parametrically  forced  uncoupled  pendulum.  A 
modal  expansion  of  the  displacement  field  reduces  the  coupled 
system  to  the  following  set  of  N  coupled  oscillators.  The  phase 
space  of  our  truncated  system  for  analysis  is  given  by  the 
2^  +  3  dimensional  vector  (T,  Z)  where 

f  =  lz.vC4})  -  {x2A'-l’X2A'})and  T  =  (6.0,  \|<3) 

Here  Z  represents  the  scaled  states  of  the  first  N  rod  modal 
oscillators,  and  m'3  represents  the  angle  of  die  periodic  forcing. 


The  2  to  1  internal  resonance  in  the  system  results  in  a  stable 
period  2  motion  in  the  pendulum  for  a  wide  range  of  parameter 
values,  including  a  very  small.  (The  linear  coupled  modes 
ex  lbit  period  1  behavior  when  the  pendulum  is  in  a  period  2 
eye  e.)  Notice  that  there  appears  to  be  no  period  doubling  route 
to  chaos  as  the  parameter  is  increased.  That  is,  the  period  2 
cycle  appears  to  undergo  instability  directly  to  a  chaotic 
attractor.  Also,  there  is  a  region  of  hysteresis  between  the  period 
2  and  the  chaotic  attractor. 

When  viewing  the  dynamics  of  extended  behavior  in  spatially 
extended  systems,  in  general  it  is  not  sufficient  to  represent  the 
dynamics  by  sampling  data  at  individual  points,  as  in  the 
method  of  time  delay  embedding  (Schwartz  and  TTiandaf, 
1996b).  instead,  the  method  of  snapshots  (Sirovich.1987)  is  a 
version  of  the  Karhunen-Loeve  method  to  analyze  the  data,  and 
computes  the  topological  dimension  of  the  attractor. 
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Figure  3.  KL  dimension  as  a  function  of  a  exhibiting 
hysteresis.  Notice  that  when  chaos  appears  in  resonance,  it  is 
high  dimensional. 


Figure  3  shows  the  dimension  as  a  function  of  the  forcing 
amplitude.  Plotted  are  two  runs  for  both  increasing  and 
decreasing  a .  For  low  parameter  values,  the  dimension  is  just  4. 
As  the  parameter  is  increased  past  the  critical  onset  of  chaos,  we 
see  a  dramatic  change  in  the  size  of  the  dimension,  which  is  now 
more  than  three  times  that  of  the  slaved  periodic  motion.  The 
system  must  have  excited  more  modes  in  order  to  generate  such 
a  high  dimension. 

Non-resonance 

If  the  system  is  tuned  away  from  resonance,  the  situation 
becomes  more  complicated.  Whereas  resonance  had  periodic 
oscillations  residing  on  a  low  dimensional  manifold,  chaotic 
oscillations  now  appear  for  both  small  and  large  g ,  which 
designates  a  weak  and  strong  coupling,  respectively.  In  Fig.  4, 
projection  of  the  motion  of  the  pendulum  onto  the  (0, 6)  plane 
shows  die  phase  portrait  when  die  dynamics  is  sampled  at  die 
forcing  period.  Notice  the  structure  that  appears  due  to  the 
manifolds  which  bound  die  motion  on  a  two-dimensional 
surface.  In  contrast,  Fig.  5  shows  die  same  projection  when  the 
coupling  between  the  rod  and  die  pendulum  is  large.  Here  the 
motion  is  not  constrained  by  any  topological  structure  since  the 
dynamics  is  now  higher  dimensional. 


Strong  Coupled  Chaos 

10  , 
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-I  •  ...  I 

°12  3  4 

f  (mod(jc)) 

Figure  5.  High  dimensional  chaos,  when  p.  =  0.288 

Examination  of  the  KL  dimension  reveals  that  although  it  is 
low  for  the  weakly  coupled  case,  it  is  not  necessarily  as  low  as 
the  resonant  case.  This  is  due  in  part  to  intermittent  escape  off 
the  manifolds.  To  see  this,  we  explicidy  compute  the  dynamics 
using  Eqs.  1  -2,  and  compare  it  to  the  dynamics  approximated  on 
the  manifold  using  a  singular  (or  center)  manifold  expansion 
(Carr,  1981),  (Schwartz  and  Georgiou,  1999).  The  error 
between  the  two  is  then  explicidy  computed  at  the  rod  tip.  For 
small  coupling,  the  error  is  small  for  most  of  the  time,  but  is  not 
uniformly  small,  since  the  scale  is  absolute.  The  dynamics  does 
in  fact  leave  the  manifold  intermittendy.  Therefore,  on  average, 
the  dynamics  remains  low  dimensional.  On  the  other  hand,  the 
strongly  coupled  case  shows  that  the  dynamics  rarely  visits  the 
low  dimensional  manifold. 
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Figure  6.  The  error  between  the  actual  and  predicted 
dynamics  of  the  rod  at  the  tip  for  weak  coupling  (left)  and 
strong  coupling  (right). 
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Figure  4.  Low  dimensional  chaos,  when  p  =  0.025  . 


EXTREME  PARAMETRIC  SENSITIVITY 

The  previous  discussion  should  make  it  clear  that  for  either 
the  near-resonant  case  or  non  resonant  case,  there  exist  multiple 
attractors.  When  the  attractors  are  chaotic  and  non-identical, 
one  needs  to  classify  the  attractors.  One  way  to  achieve  this  is  to 
compute  the  number  of  positive  Lyapunov  exponents  for  each 
attractor.  This  gives  a  statistical  measure  of  the  average  number 
of  unstable  directions  along  the  attractor.  We  have  done  this  for 
the  near-resonant  case,  and  have  found  the  existence  of  three 
distinct  attractors:  periodic,  chaotic  with  one  positive  Lyapunov 
exponent,  and  one  chaotic  attractor  with  two  positive  Lyapunov 
exponents.  (See  (Schwartz,  Wood,  Georgiou,  1999)  for 
computational  details. l  Figure  7  illustrates  the  complex 
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intertwining  in  parameter  space  of  the  three  attractors.  Such 
complexity  forces  one  to  make  large  errors  in  the  long  time 
prediction  of  the  type  of  behavior  one  expects.  That  is,  for 
arbitrary  small  amounts  of  uncertainty  in  the  parameter  value,  it 
is  almost  impossible  to  predict  the  outcome  observed  even 
statistically. 


#  of  Po*ftiv»  Lyapunov  Exponents 

Figure  7.  Number  of  positive  Lyapunov  exponents  as  a 
function  of  the  forcing  amplitude  in  the  near  resonance 
parameters. 

We  perform  the  following  uncertainty  experiment.  Pick 
parameter  a0 ,  and  measure  number  of  positive  Lyapunov 
Exponents,  NLE(aQ) .  Perturb  the  parameter  by  a  small 
amount,  e ,  called  the  uncertainty,  which  yields  =  a0  +  e .  If 

Nle( ai )  *  NLE^a o)  ’  ^en  ao  defined  to  be  uncertain .  We  then 

plot  log  of  uncertain  points  versus  the  log  of  the  uncertainty. 
The  result  is  shown  in  Fig.  8.  he  slope  is  computed  to  be  about 
0.012,  The  dashed  line  designates  what  a  good  experiment 
would  yield. 
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Figure  8.  Fraction  of  uncertainty  parameter  values  as  a 
function  of  uncertainty.  The  dashed  line  illustrates  an  ideal 
experiment  for  prediction.  In  extreme  parametric 
uncertainty,  large  increases  in  precision  (decreases  in 
uncertainty)  do  not  significantly  decrease  the  number  of 
uncertain  points.  Hence,  the  obstruction  to  predictability. 


DISCUSSION 

In  the  above  sections,  we  briefly  outlined  the  dynamics  of  a 
simple  continuum  mechanics  model  in  near-resonance  and  non- 
esonances  cases.  In  the  resonant  case,  where  the  ffequenev 
radon  between  the  pendulum  and  rod  was  1:2,  we  found  the  co¬ 
existence  between  periodic  motion  constrained  to  a  surface,  and 
high  dimensional  (m  a  KL  measure)  of  chaotic  motion.  In  the 
non-resonance  case,  we  found  low  dimensional  and  high 
dimensional  chaos  co-existing  in  large  regions  of  forcing 
amplitude.  As  a  result  of  multiple  attractors  in  the  resonance 
cases  we  defined  and  applied  a  measure  of  extreme  parametric 
sensitivity,  which  ultimately  obstructs  asymptotic  predictability 
in  certain  parameter  regions  of  interest.  . 

In  addition  to  asymptotic  attractors,  there  do  exist  temporary 
transient  phenomena  in  this  model  as  well  as  many  other 
nonlinear  systems,  but  might  be  long  lived  (Schwartz  and  ‘ 
nandaf,  1996a),  (Triandafand  Schwartz,  2000).  Such  long 
hved  (ransients  have  all  the  characteristics  of  chaos,  but  they  are 
not  attracting.  In  Fig.  9,  we  see  that  for  a  truncated  one  mode 
model  of  Eqs.  1-2,  a  chaotic  attractor  exists  at  the  right  (open 
squares).  However,  for  lower  values  of  forcing,  the  attractor 
asymptotes  to  a  period  2  attractor. 

In  many  applications,  it  is  desirable  in  the  resonance  cases  to 
keqi  the  asymptotic  periodic  state  from  occurring.  Solving  this 
difficult  control  problem  is  indeed  possible  if  one  implements  a 
topological  procedure  known  as  segmentation  (Triandafand 
Schwartz,  2000).  That  is,  the  goal  of  the  control  dynamics  is  to 
sustain  chaos  were  there  is  none  by  using  perturbations  of  die 
parameter  a .  Figure  9  illustrates  the  results  of  die  procedure, 
where  chaos  (dots)  is  sustained  over  the  whole  region  shown  In 
a  measure  theorede  sense  based  on  the  above  definition  of 
extreme  parametric  sensitivity,  the  technique  increases  ones 
predictability  by  eliminating  the  periodic  attractors. 


0.250  C-275  0.30C  0.325  0.350  0.375 


a 

I  igure  9.  Transient  and  sustained  chaos  of  a  truncated 
model 

Die  actual  parameter  fluctuations  used  to  sustain  chaos 
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throughout  the  parameter  region  are  shown  as  a  function  of 
time.  Notice  that  as  we  move  further  away  from  the  original 
chaotic  attractor,  the  frequency  of  intervention  increases. 
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Naval  Research. 
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